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A sensor was developed to quantitatively measure perturbations which change the volume of a 
wave chaotic cavity while leaving its shape intact. The sensors work in the time domain by using 
either scattering fidelity of the transmitted signals or the classical analog of the Loschmidt echo. 
The sensors were tested experimentally by inducing volume changing perturbations to a one cubic 
meter pseudo-integrable, real-world cavity. Perturbations which caused a volume change that is as 
small as 54 parts in a million were quantitatively measured. These results were obtained by using 
electromagnetic waves with a wavelength of about 5cm, therefore, the sensor is sensitive to extreme 
sub- wavelength changes of the boundaries of a cavity. The experimental results were compared 
with Finite Difference Time Domain (FDTD) simulation results, and good agreement was found. 
Furthermore, the sensor was tested using a frequency domain approach on a numerical model of 
the star graph, which is a representative wave chaotic system. These results open up interesting 
applications such as: monitoring the spatial uniformity of the temperature of a homogeneous cavity 
during heating up / cooling down procedures, verifying the uniform displacement of a fluid inside a 
wave chaotic cavity by another fluid, etc. 



I. INTRODUCTION 



Many sensors, such as SONAR, rely on direct line 
of sight information. In that case, monitoring a cavity 
which has an irregular geometric shape, including hidden 
regions, may require installing multiple sensors through- 
out the cavity for a comprehensive coverage. However, 
most real world cavities have irregular shapes. This irreg- 
ularity has the benefit of facilitating the creation of ray 
chaotic trajectories. The study of waves propagating in- 
side these ray chaotic cavities, in the semi-classical limit, 
is called wave chaos [TJ [2] . Wave chaos is essentially the 
manifestation of the underlying ray chaos on the proper- 
ties of the waves whose wavelength is much smaller than 
the typical dimensions of the cavity. For instance, ran- 
dom matrix theory has been shown to describe the spec- 
tral statistics of quantum systems with chaotic classical 
counterparts [3-6 . Ray chaos is characterized by sensi- 
tive dependence of ray trajectories to initial conditions. 
The effect of perturbations on waves propagating in such 
cavities was studied using the concept of the scattering 
fidelity [7] . Scattering fidelity is a normalized correlation 
between two cavity response signals as a function of time; 
the response signals are typically collected before and af- 
ter a perturbation to the cavity. The scattering fidelity 
decay resulting from global (as opposed to local) pertur- 
bations which change all the boundary conditions of wave 
chaotic cavities has been studied [8 -10 . Global pertur- 
bations that change only one of the walls of the cavity 
have also been considered [TTJ [12] . The scattering fidelity 
decay associated with global perturbations is either an 
exponential or Gaussian function of time. On the other 
hand, it has been shown that the scattering fidelity as- 



sociated with a local perturbation has a slower algebraic 
decay [13] . The fidelity decay induced by perturbing the 
cavity coupling has also been considered [14 . On the 
other hand, the effect of a local perturbation on a quan- 
tum mechanical system is studied using the Loschmidt 
echo concept [T5], 

Despite the extensive research on scattering fidelity, 
a practical sensing application of the scattering fidelity 
concept has not been explored. In previous work, wave 
chaotic sensing techniques that allow a comprehensive 
spatial coverage using a single sensor were introduced 
[T6HT8] . These techniques rely on the wave chaotic nature 
of most real world cavities. When a wave is broadcast 
into a cavity to probe it, the response signal consists 
of reflections that bounce from almost all parts of the 
cavity; this is due to the underlying spatial ergodicity of 
ray trajectories in ray chaotic cavities. Therefore, the 
response signal essentially "fingerprints" the cavity, and 
it enables the detection of changes to the cavity. 

The wave chaotic sensing techniques developed in 
Refs. pj)HT8] were not used to quantify any kind of per- 
turbation. Local and global perturbations to the bound- 
aries of the cavity, and perturbations to the medium of 
wave propagation within the cavity, were all shown to be 
detectable [17]. However, the quantification of a pertur- 
bation was not accomplished. On the other hand, a re- 
markably sensitive quantification of a perturbation which 
involved translation of a sub- wavelength object over sub- 
wavelength distances was successfully demonstrated [19] . 
However, the quantification was based on an empirical 
law that is specific to the system and perturbation at 
hand. This is because the effect of the perturbation on 
the dynamics of the waves propagating inside the wave 
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chaotic cavity is not straightforward [19 . In this pa- 
per, we focus on a single class of perturbation whose ef- 
fect can be theoretically predicted, and we propose two 
time domain techniques to measure that particular kind 
of perturbation in any cavity. 

In this paper, we focus on quantifying volume changing 
perturbations (VCP) to a wave chaotic scattering sys- 
tem. A VCP changes the volume of a cavity, but it 
may slightly change its shape as well. A special kind 
of VCP is a volume changing and shape preserving per- 
turbation (VCSPP). In Sec.[nJ the theoretical prediction 
of the effect of VCSPPs is discussed. Sec. [TT] proposes 
two time domain techniques to quantify VCPs. As in 
previous work, these techniques are based on the scat- 
tering fidelity and the classical analog of the Loschmidt 
echo [T6HT8] . Sec. Ill presents the experimental test of 
these two VCSPP sensing techniques, along with a head- 
to- head numerical validation. The experimental test is 
carried out inside a pseudo-integrable, real world cavity 
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using electromagnetic waves. Sec. [TV] provides a test of 
the sensing techniques in a numerical model of the star 
graph, which is a quasi-lD wave chaotic system. Sec. [IV 
also shows the relative merits of approaching the prob- 
lem in the frequency domain. Sec. |V| discusses practical 
applications of the VCSPP sensor, and Sec. VI provides 
a conclusion. 



II. THEORY AND APPROACH 

Consider a generic wave chaotic cavity with volume V\ , 
which is considered as a baseline system (see Fig. [TJa)). 
The schematics in Fig. [TJa)&(b) illustrate the cavity as 
a stadium billiard, but the cavity is considered generic 
throughout Sec. [TTJ Suppose that the baseline cavity is 
perturbed such that each of its three length dimensions 
increase by a factor of P. This amounts to a VCSPP, by 
a factor of P 3 ; the perturbed cavity has a volume of V2 = 
P 3 l/i (see Fig. [IJb)). Fig. [jja)&(b) show a brief pulse 
being broadcast into the cavity. The response signal to 
the pulse is called the sona. The typical durations of 
the pulse and the sona are Ins and 10/is respectively. 
The sona from the baseline cavity (which is referred to 
as baseline sona) and the sona from the perturbed cavity 
(which is referred to as perturbed sona) are expected to 
be related, under certain conditions which are discussed 
in Sec. |III E| For instance, if P > 1, a signal feature in 
the perturbed sona is expected to be delayed by a factor 
of P compared to its appearance in the baseline sona. 
The sensor is designed to enable the measurement of the 
value of P, which effectively quantifies the VCSPP, by 
using the theoretically predicted effects of the VCSPP 
on the dynamics of the waves. Another practically useful 
capability of the sensor is to check if the perturbation is 
indeed a VCSPP, and not just merely a VCP. 

Next, consider the scattering parameters of the cavi- 
ties as a function of frequency. The |6i2| 2 as a function 
of frequency of the baseline and the perturbed 2-port 




FIG. 1: Schematic illustrating Volume Changing & Shape 
Preserving Perturbations (VCSPP). (a) Sona is collected from 
a baseline cavity of volume V\. (b) Sona is collected from a 
perturbed cavity of volume V2. (c) Pulse exciting the reso- 
nances of the baseline system, (d) The same pulse exciting 
the perturbed resonances of the perturbed system. 



cavity are schematically shown in Figs. []Jc)&(d). We 
expect a precise mathematical relationship between the 
scattering parameters of the baseline and perturbed cav- 
ities as a function of frequency. In particular, if P > 1, 
the baseline spectrum can be obtained by stretching out 
the perturbed spectrum by a factor of P along its fre- 
quency axis. This is precisely the prediction about the 
effect of VCSPPs on the dynamics of the waves. This pre- 
diction can be used to measure the perturbation in the 
frequency domain. As opposed to the resource intensive 
frequency domain sensing, it is practically preferable to 
use a time domain interrogation of the baseline and the 
perturbed cavity by measuring the sonas. However, it is 
useful to look at the problem in the frequency domain to 
understand the limitations of the time domain approach 
discussed in Sec. MI El 

There are two classes of time domain sensing tech- 
niques that can be used to quantify VCSPPs. The first 
technique relies on the scattering fidelity [7]. Consider 
two sonas X and Y, which are real voltage versus time 
signals with zero mean values. The scattering fidelity 
(SF) of X and Y is simply their Pearson's correlation as 
a function of time, t [7]; 



SF(t) = 



EZ=t +At X[m]Y[m] 

•=* +At y[ m ]2 



Z-^m—t L J 

lYZl\ +At x[m?YZl\ 



(1) 



where At is typically chosen to be the time it takes the 
waves to traverse the cavity, at the very least, once (i.e. 
in order of magnitude of the ballistic flight time). The 
SF(t) of two sonas can take real values ranging from 1 
(i.e. perfect correlation at time t) to —1 (i.e. perfect 
anti-correlation at time t). If SF(t) is 0, then the sonas 
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are not correlated at time t. The SF of the baseline 
and the perturbed sonas is not expected to stay close to 
1 throughout time. However, the SF of the perturbed 
sona and the baseline sona whose time axis is scaled us- 
ing the optimum stretching/squeezing factor is expected 
to stay close to 1 throughout time. The optimum stretch- 
ing/squeezing factor is expected to be equal to P, which 
is also related to the magnitude of the perturbation. 

The second technique to quantify VCSPPs utilizes the 
classical analog of the Loschmidt echo through classical 
time reversal mirrors [20] . To see the operation of a time 
reversal mirror, consider a two port cavity. Suppose that 
a pulse is broadcast into the baseline cavity through port 
1, and a baseline sona is recorded through port 2. If the 
baseline sona is time reversed and broadcast back into 
the cavity through port 2, a time-reversed version of the 
original pulse reconstructs at port 1. The reconstructed 
pulse approximates the time reversed version of the orig- 
inal pulse broadcast into the cavity. However, if the time 
reversed baseline sona is broadcast into a perturbed cav- 
ity, then the reconstructed pulse will more poorly approx- 
imate the time reversed version of the original pulse. The 
time axis of the baseline sona needs to be scaled using the 
optimum factor before it is time reversed and broadcast 
into the perturbed cavity; this is assuming that the per- 
turbation is VCSPP. The optimum stretching/squeezing 
factor is expected to result in a reconstructed pulse that 
best approximates the original pulse. Once again, the 
optimum stretching/squeezing factor is expected to be 
P. 

Time-reversal mirrors have found a wide range of prac- 
tical applications such as crack imaging in solids [2TJ [22] , 
and improved acoustic communication in air [23] , among 
other things [20j [24]. Time-reversal mirrors have been 
shown to benefit from the cavity's underlying ray chaos, 
which is prevalent in most real world cavities [25]. Re- 
cently, it was proposed that time reversal mirrors could 
also be applied to quantum systems [26] . The robustness 
of time-reversal mirrors in a scattering medium undergo- 
ing perturbations has also been studied [27] . 



III. TESTING THE QUANTITATIVE SENSOR 
IN THE TIME DOMAIN 



A. Experimental Setup 

The cavity that is used to test the sensing techniques is 
an approximately lm 3 (i.e. dimensions of 1.27m x 1.27m 
x 0.65m) aluminum box that has scatterers and interior 
surface irregularities which facilitate the creation of ray 
chaotic trajectories (see Fig. [2|. The cavity is pseudo- 
integrable because it has parallel walls which may sup- 
port integrable modes in addition to the chaotic modes. 
Overall, the cavity represents a real world case in which 
the sensor would operate. There are two ports that con- 
nect the cavity to a microwave source and an oscillo- 
scope. Each port consists of a monopole antenna of 



length w 1cm, and diameter « 1mm. The monopole 
antennas are mounted on two different walls of the cav- 
ity. An electromagnetic pulse with a center frequency of 
1GHz, and a Gaussian envelope of standard deviation 
Ins is typically broadcast into the cavity through port 
1. The resulting sona signal is collected at port 2 by the 
oscilloscope, and it is digitally filtered to minimize noise. 

Experimentally inducing a VCSPP can be more chal- 
lenging than inducing a VCP, which may slightly change 
the shape of the enclosure. VCSPPs can be realized by 
changing the speed of wave propagation within the cavity 
(i.e. changing the electrical volume). The speed of light 
(for the electromagnetic experiment at 7GHz) within the 
cavity can be changed by filling up the cavity with differ- 
ent gases which have similar dissipation and dispersion 
properties. For example, the relative dielectric constant 
(e r ) of air (at « 50% relative humidity), nitrogen gas, 
and helium gas are e r ^ a i r = 1.000635, e r ,N 2 — 1.000547, 
e r ,He = 1.000065, respectively, at a temperature of 20°C 
and a pressure of 760torr [28j [29] . As discussed later in 
Sec. |IIIB[ it was shown that the slightly different dis- 
sipation values of these gases does not change the sona 
signal in any perceivable way. However, the slightly dif- 
ferent speed of light values of these gases was seen to 
significantly change the sonas collected. 

The experimental procedure for filling up the cavity 
with different gases is as follows (see Fig. |2|. There is 
a gas inlet on the top wall of the cavity. The gas inlet 
was connected to a gas tank via a plastic tube and a 
long copper tube coil. The long copper tube coil allowed 
the gas to reach room temperature before it gets into the 
cavity. There was a pressure regulator in between the gas 
tank and the copper tube to control the rate of flow of 
the gas. There were three gas outlets near the top wall 
of the cavity, and three gas outlets near the bottom wall 
of the cavity. The diameter of the gas inlet and outlets 
was about a fifth of the wavelength, so that no significant 
microwave leakage occurred. Depending on the density 
of the gas which was being pumped into the cavity, half of 
the outlets (near the top or bottom wall) were closed off 
with tape. This procedure helped to displace the existing 
gas and retain the gas being pumped into the cavity. 

As detailed in Sec. IIID[ sona signals were collected 
from the cavity both during and after the gas transfer 
process. The sona signals that were collected during the 
gas transfer indicate when the cavity is almost fully filled 
with the new gas. The sona signals that were collected 
after the gas transfer were used to measure the VCSPP 
that was induced. 



Finite Difference Time Domain (FDTD) 
Simulations 



The experiment described in Sec. |III A| was modeled 
using a Finite Difference Time Domain (FDTD) code. 
The FDTD code solves Maxwell's Equations inside a 3D 
numerical model of the experimental cavity introduced 
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FIG. 2: Schematic of the experimental set up to induce VC- 
SPPs in an electromagnetic cavity, and the equipment needed 
to implement an electromagnetic time reversal mirror. The 
VCSPP is induced by filling the cavity with helium or nitro- 
gen gas. The gas transfer is carried out using a long copper 
tube that helps to warm the gases up to room temperature. 
There is a gas inlet, and six gas outlets on the walls of the 
cavity. The cavity has two antennas that are connected to 
a microwave source and an oscilloscope. The electromagnetic 
time reversal is carried out as follows. First, the original pulse 
is broadcast through antenna 1 (1), and the resulting sona 
is collected at antenna 2 (2). Next, the time reversed sona 
is injected into the system at antenna 1 (3) utilizing spatial 
reciprocity to retrieve the reconstructed time reversed pulse 
at antenna 2 (4). Experimental data are shown for each step. 



in Sec. |III A| The code for this simulation was optimized 
for parallel computers and for the simulation of rever- 
beration chambers [30] [31] . The FDTD simulation of the 
cavity enabled a direct comparison of experimental and 
simulation results. 

The FDTD simulation formed a 3D model of the ex- 
perimental cavity by using spatial cubic cells with an 
edge length of Ax w 3.49mm, and the cavity consisted 
of 364 * 364 * 188 cells. The smallest time step taken to 
propagate solutions of Maxwell's Equations through the 
cells was At = 6ps. Therefore, the Courant's number for 
computational stability was C J^^ X ~ 0.89 < 1 for all the 
media (e r ) considered, where c is the speed of light in vac- 
uum [32 . The model of the cavity also had two ports, 
with antennas that are similar to the monopoles used 
in the experiment. The electromagnetic pulse broadcast 
into the model was also similar to the one broadcast ex- 
perimentally (i.e. center frequency of 1GHz, and width 
of Ins). The maximum and minimum of the ratio of the 
wavelength to the cubic cell dimension were 14 and 11 
respectively. 

In the experiment, the main source of dissipation is 
ohmic loss from the aluminum walls of the cavity. In 
the simulation, the walls were assumed to be lossless for 



simplicity. Instead, an equivalent loss was introduced 
within the medium of wave propagation to achieve the 
same quality factor as the experimental cavity [30 . To 
accomplish this, a uniform conductivity of 10 _5 5'/m was 
introduced throughout the interior of the cavity model. 

As mentioned in Sec. Ill A[ the experiment relies on 
changing the electrical volume (at 7GHz) within the cav- 
ity by changing the gases filling the cavity. The gases 
used were air, nitrogen gas, and helium gas. For typical 
laboratory atmospheric conditions, the electromagnetic 
loss in air at 1GHz can be mainly attributed to oxygen 
molecules [30] [33]. The specific attenuation of oxygen 
is 0.007gLB I km whereas the specific attenuation of water 
vapor is just 0.003dB/km at 7GHz. The conductivity of 
air is estimated to be 4.3 * 10 _9 5'/m [30]. Therefore, the 
dissipation inside the cavity filled by helium gas or nitro- 
gen gas was roughly modeled by introducing an equiva- 
lent loss of 10- 5 S7m-4.3*10- 9 5ym = 9.9957* 10~ 6 S/m 
uniformly throughout the cavity (once again, the walls 
were assumed to be lossless in the model). It was shown 
that sonas that are collected from the cavity model with 
10 _5 5'/m conductivity, and sonas that are collected from 
the cavity model with 9.9957* 10 _6 S'/m are identical (i.e. 
they have a fidelity of 1). This simulation result proved 
that the difference in loss among air (at « 50% relative 
humidity), nitrogen gas, and helium gas is not significant. 

As will be discussed later in Sec. |IIID[ the gases that 
experimentally filled the cavity are not pure. The effec- 
tive e r of the gases that filled the cavity are different from 
the e r values for pure gases. The simulation of the cav- 
ity that is filled with helium, nitrogen or air is done by 
using the effective e r values. This allows the simulation 
to better model the experimental reality. The e r values 
used to simulate a cavity filled with air, nitrogen and 
helium are 1.000576, 1.000547, and 1.00017 respectively. 
These values were chosen to match the experimental re- 
sults that will be discussed in Sec. MIDI The values were 
determined by anchoring the effective e r of nitrogen to 
the value for pure nitrogen gas, and finding the effective 
e r values for the other gases to match the experiment. 



C. Sensitivity of the Quantitative Sensor 

Here, the sensitivity of the sensor is described by deriv- 
ing an expression for the minimum volume changing per- 
turbation that can be measured. Then, it will be shown 
that the VCSPP which is induced when either of the 
three gases in the experiment (i.e. air, nitrogen and he- 
lium gas) is displaced by another one, can be measured 
using our experimental set up. 

Consider using the sensing technique based on scat- 
tering fidelity, which was introduced in Section [TTJ The 
technique relies on comparing the baseline and perturbed 



sonas as a function of time. In Sec. IIII CI it is assumed 



that the perturbation increases the volume by a factor 
of P 3 , where P > 1. At time t = 0, the sonas are ex- 
pected to be similar for small enough At (see Eq. [j), 



hence SF(t = 0) « 1. At any other time t, there may be 
a perceptible difference between the sonas. Any particu- 
lar signal feature in the baseline sona at time t is expected 
to be seen in the perturbed sona at time t+t gap (t), where 
tgap(t) = P t — t; here, t gap (t) is defined as the time gap 
that develops between two identical features in the base- 
line and perturbed sona at time t, where t is measured 
within the baseline sona. This is because the baseline 
sona stretched out along its time axis by a factor of P 
should approximate the perturbed sona, as discussed in 
Sec. ED 

The minimum volume changing perturbation that can 
be quantified is the minimum of P 3 = ( t+tg t ap( ^ ) 3 . On 
the other hand, when the SF of the baseline and per- 
turbed sona is computed at any time t (see Eq. [T]) , there 
are two necessary conditions that should be satisfied in 
order to be able to measure the perturbation. First, the 
signal-to-noise-ratio (SNR) of the baseline sona at time 
t should be well above 1. This is because the SF of the 
baseline and perturbed sonas (when the SNR is close to 
1) would simply be the correlation of two noisy signals. 
Second, t gap (t) should be, conservatively, greater than 
half of the period of the oscillations in the sona signals. 
Otherwise, if t gap (t) is much smaller than half a period of 
the sona oscillations, the SF(t) will not be convincingly 
lower than its maximum value at t = (which is ~ 1 
for the appropriate At value in Eq. [T]) , hindering reliable 
measurement of the perturbation. These two conditions 
guarantee that the SF(t) of the baseline and perturbed 
sona can be used to measure the perturbation. The mini- 
mum of P 3 = ( t+tfl ° p ^ ) 3 subject to these two conditions 
is (1 + Tf, — — ) 3 , where T is the period of sona oscil- 

V J- NoiseLevel ' 1 

lations and T^oiseLevei is the time at which the SNR of 
the sona approaches 1. 

Therefore, the minimum perturbation that can be 

quantified by the sensor (i.e. (1 + T T/2 ) 3 ) depends 

on the wavelength of the waves used to prooe the cavity, 
the dissipation in the cavity and the SNR of the system. 
The SNR in turn depends on the noise in the system, 
and the dynamic range of the wave generation and detec- 
tion equipment. The equipment used in this experiment 
are shown in Fig. [2] For the electromagnetic experimen- 
tal set up that is probing the Ira 3 cavity discussed in this 
work, using 5cra wavelength waves, a VCSPP that is as 
small as 4 parts in 10 5 can be measured. 

Getting back to the experimental set up shown in 
Fig. |2j the VCSPP was induced by changing the gas fill- 
ing the cavity from air to nitrogen gas or to helium gas. 
This results in « 0.05% change in the e r of the gases, 
which is equivalent to a VCSPP of at least 13 parts in 
10 5 . Therefore, the experimental system is expected to 
detect the change in electrical volume induced when one 
of these gases displaces the other inside the cavity (as- 
suming that the gases are pure). 




Time (jis) 



FIG. 3: Scattering fidelity of sona signal from a cavity that 
is filled with air (baseline) , and sona from a perturbed cavity 
that had 2, 4, 6, 8, and 10 minutes of nitrogen gas inflow. The 
perturbation gets stronger as the concentration of nitrogen 
gas increases in the perturbed cavity. 



D. Results from the Experiment and the FDTD 
Simulation 

1. Sensing Using Scattering Fidelity 

The air in the cavity at ~ 20° C and « 50% relative 
humidity, was systematically displaced with nitrogen gas 
at room temperature. The nitrogen gas was pumped into 
the cavity at 30psi gauge pressure as the air flowed out 
through the gas outlets of the cavity. Every two minutes, 
the flow was stopped, and 10 nominally identical sonas 
(which are actually almost identical) were measured from 
the cavity, and these were averaged together. The aver- 
aging is done after aligning the sonas to eliminate the 
adverse effects of trigger jitter in the data acquisition 
system. In this manner, five averaged sonas were col- 
lected from the cavity as the cavity was filled with more 
and more pure nitrogen gas. Each of these five sonas 
were compared with a sona that was collected from the 
original cavity filled with air. The comparison was done 
by computing the scattering fidelity (see Eq. [I]) of the 
sona from airy cavity and a sona from a partially air- 
filled cavity. Fig. [3] shows these scattering fidelities. The 
concentration of nitrogen increases with the number of 
minutes of nitrogen inflow into the cavity. Therefore, 
Fig. [3] shows scattering fidelities of VCSPPs which get 
progressively stronger. 

The scattering fidelity of a VCSPP shows oscillation 
whose period is inversely related to the strength of the 
perturbation (see Fig. |3|. The oscillation in the scatter- 
ing fidelity can be explained by specifically examining the 
scattering fidelity of sona from the air filled cavity, and 
sona from the nitrogen gas filled cavity (see Fig. [4^a)). 
Here, both of the sonas were obtained by averaging over 
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FIG. 4: (a) Examining the scattering fidelity oscillation for 
VCSPPs using sona from air filled cavity and sona from ni- 
trogen filled cavity. Each of the sonas are averaged over 100 
sona samples, (b) The sonas near t = have fidelity of 1. (c) 
The sonas are out of phase by half a period around t = 4/xs. 
(d) The sonas are out of phase by a period around t = 8 /us. 
Besides, the phase shift between the envelopes of the sonas 
becomes significant. 



100 nominally identical sona samples. For times rela- 
tively close to t = 0, the difference in the speed of light 
between air and nitrogen does not show up when signals 
that traveled through these two gases are compared (see 
Fig. [4^b)). However, at t « 4/is, a phase shift of half 
a period (i.e. 0.07ns for the probing pulse centered at 
1GHz) develops between the two signals (see Fig. |4^c)). 
Thus, the scattering fidelity between the two sonas be- 
comes ~ —1. The relative phase shift between the two 
sonas increases to one full period at about 8/is, and hence 
the scattering fidelity recovers to almost 1 (see Fig.[2Jd)). 
However, as can be seen in Fig. [4^d), in addition to the 
phase shift that develops between the individual oscil- 
lations of the two sonas, a relative phase shift starts to 
develop between the envelopes of the sona signals. Thus, 
the scattering fidelity of VCSPPs is not expected to oscil- 
late between 1 and —1 indefinitely. Rather, it is expected 
to decay while oscillating [9]. However, we do not an- 
ticipate seeing this fidelity decay with our measurement 
system for this particular perturbation because the SNR 
of the sonas approaches unity after about 10/is, and the 
fidelity does not decay appreciably within this time for 
this particular perturbation. 

The speed of light at 1GHz in pure nitrogen gas was 
faster than it was in the laboratory air (which had about 
50% relative humidity). Therefore, as discussed in Sec.|H| 
the VCSPP can be quantified by finding the optimum 
stretching factor to be applied on the sona that is col- 
lected from the cavity filled with nitrogen. The goal is to 
recover the scattering fidelity of the stretched "nitrogen- 
sona" and the "air-sona" to 1. This was achieved by 
using a stretching factor given by (4/is + 0.07ns)/4/is; 



because, t gap = 0.07ns half a period of oscillation 
of the sona) at t = 4/is (SF « — 1) based on the dis- 
cussion in Sec. |III C| The resulting scattering fidelity is 
plotted in red (see Fig. [5] ) with the scattering fidelity of 
the unmodified sonas which is shown in blue. The opti- 
mum stretching factor is (4/is + 0.07ns) /4/is « 1.000018. 

Whereas, 4 / e J r,Air ~ 1.000044 is the expected value of 
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P. In other words, ideally a 44ppm (parts per million) 
change across each electrical dimensions of the cavity is 
expected, however, an ISppm change is measured. 

The discrepancy can be explained by the fact that the 
nitrogen gas that filled the cavity has an effective e r that 
is probably much closer to the effective e r of air. Fig. [3] 
shows that it can take several minutes to displace the air 
in the cavity with pure nitrogen. Therefore, we do not 
expect the e r of the nitrogen gas that filled the cavity 
to be the same as the literature value for pure nitrogen 
gas. The same conclusion applies to the other gases fill- 
ing the cavity. There are several reasons for the expected 
discrepancy between the literature value of e r and the ef- 
fective e r of the gases filling the cavity experimentally. 
The helium and nitrogen gases are of industrial quality, 
which is not perfectly pure. Besides, the cavity is not nec- 
essarily air tight; this could lead to leakage of air into the 
cavity when the cavity is not pressurized. Finally, the e r 
of the gases is a function of temperature and relative hu- 
midity (for the case of air) [29]. Therefore, deviations of 
laboratory temperature and relative humidity from 20° C 
and 50% could impact the effective e r value of the gases. 

The fact that the shape of the cavity was preserved 
(during the displacement of the air by nitrogen gas) can 
also be seen from Fig. [5] If the shape of the cavity were 
not preserved, it would not be possible to recover the 
scattering fidelity of the two sonas through simple nu- 
merical stretching of one of the sonas, along the time 
axis. Hence, displacing air with nitrogen gas is not just 
a VCP, but is also a VCSPP. Later, a VCP induced by 
displacement of air with helium, which is not VCSPP, 
will be discussed. 

The fidelity decay of VCSPPs, which is expected to 
be superimposed on the fidelity oscillations, can be seen 
for a stronger VCSPP. The VCSPP should be strong 
enough to bring about a significant phase shift between 
the envelopes of the sonas before their SNR deteriorates. 
For our experimental set up, such a strong VCSPP can 
be achieved by displacing the air in the cavity by he- 
lium gas. The scattering fidelity of sona from a cavity 
that is filled with air and sona from a cavity that is 
filled with helium is plo tted i n Fig. [6] in blue. Based 



on the definition in Sec. IIIC 1 t gap = 0.07ns (which is 



half the period of sona oscillation) at time t = 0.35/is. 
Thus, the optimum stretching factor was chosen to be 
(0.35/is + 0.07ns)/0.35/is = 1.0002. This optimum value 
maximizes the average value of the resulting scattering 
fidelity. Since the speed of light at 7GHz is higher in 
helium than in air, it was the "helium sona" that was 
stretched out along its time axis. The scattering fidelity 
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FIG. 5: Undoing the effect of a VCSPP. The sona that was 
collected from the cavity filled with nitrogen was stretched out 
optimally to recover the scattering fidelity to 1 throughout 
the times when the SNR of the sonas is robust. The SNR 
decreases roughly by a factor of a thousand between Ofis and 
10/is. The optimum stretching factor quantified the VCSPP. 
The fact that the scattering fidelity was recovered proved that 
the perturbation was VCSPP. 



FIG. 6: Scattering fidelity of sona from a cavity filled with 
air and sona from a cavity filled with helium. The fidelity 
oscillations and decay can be seen. The effect of the VCP 
could not be completely undone by stretching the "helium 
sona", so, the perturbation is not VCSPP. The buoyant force 
of helium gas can slightly change the shape by flexing the 
walls of the cavity. 



of the stretched "helium sona" and the "air sona" is 
also shown in Fig. [6] in red. Once again, the stretch- 
ing factor (0.35/is + 0.07ns)/0.35/is = 1.0002 approxi- 

'— w 1.000285. In other words, a 285ppm 



mates 

change across each electrical dimensions of the cavity is 
expected, and a 200ppm change is measured. This shows 
that the change in electrical volume which was induced 
by replacing the air in the cavity with helium was quan- 
tified successfully, considering the fact that the effective 
e r of the helium and air gases are probably closer than 
expected. 

However, unlike the case in Fig. [5j the effect of the 
VCP could not be undone perfectly. The scattering fi- 
delity of the stretched "helium sona" and the "air sona" 
was not close to 1 throughout time; instead it shows a 
fidelity decay. The scattering fidelity of "helium sona" 
and "air sona" is expected to oscillate between 1 and —1 
(as the phase shift between the fast oscillations of the 
sonas increases in a similar fashion to the illustration in 
Fig. Qb-d)) and decay to (as the phase shift between 
the complicated envelopes of the sonas increases). This 
decay is seen experimentally in Fig. [6j However, there 
must be another fidelity decay superimposed on the fi- 
delity decay that can be attributed to a VCP; because 
the fidelity decay could not be undone by numerically 
stretching out one of the sonas. 

A FDTD simulation of the experiment was performed 
to better understand the form of the scattering fidelity 
when comparing "helium sona" and "air sona". As dis- 
cussed in Sec. IIIB[ differences in the dissipation of he- 
lium gas and air are so minute that they do not need to be 



considered. The scattering fidelity of "helium sona" and 
"air sona", which were obtained from the FDTD model 
by broadcasting the 7GHz pulse used experimentally, are 
shown in Fig(7[ The simulation results show that the ef- 
fect of the VCSPP can be undone by applying the opti- 
mum stretching factor (i.e. (0.35/is + 0.07ns) /0.35/is = 
1.0002) to the "helium sona". The scattering fidelity of 
the "air sona" and the stretched "helium sona" is shown 
in red in Fig(7| it is close to 1 which shows that the effect 
of the perturbation can be undone by simple numerical 
stretching. The fidelity is close to 1 for times where the 
numerical errors in the FDTD are negligible. 

The difference in the results in Fig. [6] and Fig. [7| shows 
that experimentally displacing air with helium induces 
another kind of perturbation other than a volume chang- 
ing perturbation. The additional perturbation was not 
seen when nitrogen was pumped into the cavity at the 
same pressure setting (i.e. SOpsig) as was used for pump- 
ing helium into the cavity. This discrepancy can be ex- 
plained by the fact that helium gas exerts a significant 
buoyancy force which slightly flexes the walls of the cav- 
ity (about lm 2 area, and 3mm thick aluminum sheets). 
From previous work on this cavity, it was shown that a 
flexing of one of the walls of the cavity can cause a signifi- 
cant shape changing perturbation [34] . The result shown 
in Fig. [6] demonstrates that it is possible to verify if the 
shape of the cavity remained intact while its electrical 
volume changed. This verification can be simply done 
by checking if the scattering fidelity can be recovered to 
1 throughout time. The capability to detect changes to 
the shape of the cavity during a volume changing pertur- 
bation (which could be induced by a spatially uniform 
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FIG. 7: Scattering fidelity of sona from a cavity filled with 
air and sona from a cavity filled with helium. The sonas are 
generated from the FDTD model of the cavity. 
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"Squeezing" Factor 



FIG. 8: Plot of peak-to-peak amplitude of time-reversed pulse 
in a cavity filled with Helium vs "squeezing" factor used to 
scale the time axis of the sona which is collected from a cavity 
filled with Nitrogen. 



heating or cooling of a homogenous cavity) can have sev- 
eral applications as was pointed out in Sec. [TJ 

2. Sensing Using Time Reversal 

So far, only one of the sensing techniques introduced 
in Sec. [IT] is demonstrated to quantify volume changing 
perturbations. In addition to the scattering fidelity tech- 
nique, time reversal mirrors can be used to quantify vol- 
ume changing perturbations. When a sona signal that is 
collected from the cavity filled with nitrogen is time re- 
versed and broadcast into the cavity filled with helium, 
the reconstructed time reversed pulse is not expected to 



be ideal. This is because the time reversed sona would be 
traversing an effectively smaller cavity. Note that in this 
particular case, the baseline cavity is electrically larger 
than the perturbed cavity. Based on the discussions in 
Sec. [TTJ the baseline sona should be squeezed along its 
time axis using an optimum factor to recover the maxi- 
mum amplitude time reversed pulse. When the optimally 
squeezed "nitrogen sona" is time reversed and broadcast 
into the cavity filled with helium, the reconstructed time 
reversed pulse is expected to better approximate the orig- 
inal pulse. 

The improvement in the quality of the time reversed 
reconstructed pulse can be measured in various ways [17] . 
Here, the simplest measure (the peak-to-peak-amplitude 
(PPA) of the reconstructed pulse) of the quality of the 
time reversed pulse is used. Fig. [8] shows the PPA of 
the time reversed reconstructed pulse obtained when the 
"nitrogen sona" is numerically squeezed along its time 
axis by varying amounts. The optimum squeezing factor 

of 0.9998 approximates y^f^ ~ 0.99976. This means 

that a 240ppm change is expected, but a 200ppm change 
across each electrical dimensions of the cavity was mea- 
sured. Once again, the small discrepancy can be ex- 
plained by the fact that the gases in the cavity are not 
perfectly pure. This shows that the time reversal tech- 
nique can also be used to quantify volume changing per- 
turbations. 

Time reversal mirrors can also be used to detect when a 
volume changing perturbation slightly changes the shape 
of the cavity. In this case, when the "nitrogen sona" is 
time reversed and broadcast into the cavity filled with ni- 
trogen, the PPA of the reconstructed time reversed pulse 
was about 1.13V. However, when an optimally squeezed 
"nitrogen sona" is broadcast into the cavity filled with 
helium gas, the optimal value of the PPA is only about 
0.85V. This indicates that the perturbation, which is in- 
duced when the nitrogen gas inside the cavity is displaced 
with helium gas, is not just a volume changing pertur- 
bation but also a shape changing perturbation. Once 
again, this is due to the relatively strong buoyant force 
of helium gas which can flex the walls of the cavity. 

To summarize, either of the two time domain sensing 
techniques introduced in Sec. [il] can be used to identify 
and quantify a VCSPP. However, the sensing technique 
based on time reversal can have an advantage because it 
is computationally cheaper [T7] . 

E. Limitation of the Time Domain Approaches to 
Measure VCSPP 

Now we can describe the time domain approach fur- 
ther, and point out its limitations. When a cavity is 
monitored for a VCSPP, a pulse is periodically broad- 
cast into it, and a sona is collected. It is generally pre- 
ferred to use the same probing pulse (i.e. the same cen- 
ter frequency, bandwidth, amplitude, and shape) so that 
only the changes in the system show up when looking 
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at the sonas. When a VCSPP (with P > 1) occurs, 
the transmission spectrum |5i2| 2 of the system shrinks 
along the frequency axis by a factor of P as shown 
in Figs. [TJc)&(d). The probing pulse frequency cover- 
age is schematically shown in the frequency domain in 
Figs. [TJc)&(d). The baseline and the perturbed sonas 
are a result of the probing pulse exciting resonances of 
the cavity. The resonances excited by the probing pulse 
in the baseline and perturbed cavity are not all the same. 
Suppose that there is a significant overlap between the 
resonances excited by the probing pulse in the baseline 
and perturbed cavity. Under this condition, we expect 
that the baseline sona can be numerically stretched out 
by a factor of P along its time axis to approximate the 
perturbed sona. 

However, if there is no overlap between the resonances 
excited in the baseline and perturbed system, the time 
domain sensing techniques are not expected to work. 
Thus, both of the time domain sensing techniques (based 
on scattering fidelity and time reversal mirrors) face a 
limitation regarding the maximum VCSPP that can be 
measured. This limitation can be improved by increas- 
ing the bandwidth of the pulse that is used to probe the 
cavity. Doing so would effectively make the time domain 
technique closer to the frequency domain interrogation of 
the system, which does not face a limit on the maximum 
VCSPP that can be measured. In a similar spirit, using 
a pulse with a flat frequency response (such as a chirp) 
may also be helpful. 

The regime of VCSPP strengths that cannot be mea- 
sured using the time domain techniques is estimated as 
follows. Suppose that the probing pulse excites reso- 
nances with frequencies ranging from f min to / max , with 
the center frequency at / = f min +f™- ax ? and 3dB band- 
width of fmax — fmin- It is assumed that the probing 
pulse excites several resonances of the baseline system. 
In other words, f max - f min > A/, where Af is the 
mean spacing between the resonant frequencies of the 
baseline cavity. The VCSPP changes the volume of the 
cavity by a factor of P 3 , where P > 1 for simplicity. Once 
again, the VCSPP has the effect of scaling the \Si 2 \ 2 (f) 
of the cavity along the frequency axis by a factor of P. 
If ^ m p x < fmin, then the VCSPP is not expected to be 
measured by using this probing pulse in the time domain 
because there would be significant overlap between the 
resonances excited in the baseline and perturbed cavities. 

In the experiments discussed in Sec. |III A[ the prob- 
ing pulse excited resonances between f min = 6.5GHz 
and fmax — 7.5GHz. The mean spacing between the 
resonant frequencies of the baseline cavity is given by 
Af = ^rg^ ~ 22kHz; where V is the volume of the 
baseline cavity, / = 1GHz is the center frequency of the 
pulse, and A is the wavelength. Therefore, a VCSPP 
where P > ^ max ~ 1.15 is not expected to be measured 



vs the air) are both far below the maximum perturbation 
strengths that can be measured using the time domain 
techniques. 

The case of a strong VCSPP perturbation, which can- 
not be measured using the time domain techniques, is 
best considered using a simulation tool that can be eas- 
ily interrogated in the frequency domain as well. Sec. [IV 
discusses such strong perturbations, and shows how they 
can be quantified using a frequency domain approach. 



IV. SENSING USING FREQUENCY DOMAIN 
INFORMATION 



In Sec.[TTJ it was mentioned that VCSPPs can be quan- 
tified using information obtained in the time domain or 
frequency domain. The time domain approach is gen- 
erally more practical in applications. However, the fre- 
quency domain approach does not have limitations on 
the maximum perturbation value that can be quantified. 
This frequency domain approach is used to measure VC- 
SPPs on a quasi- ID system called the star graph. We 
use the star graph because it is a type of quantum graph 
that has generic properties of wave chaotic systems [35] , 
but is relatively simple to understand. It is also compu- 
tationally cheaper to simulate than the 3D wave chaotic 



frr 



in the time domain experiments discussed in Sec. Ill A| 
Clearly, the strongest VCSPP that was achieved in the 
laboratory (i.e. P ~ 1.000018 for nitrogen gas vs the air), 
and the strongest VCP (i.e. P ~ 1.0002 for helium gas 



system discussed in Sec. IIIB Besides, the star graph 
can be directly implemented in the frequency domain as 
discussed in Sec. IIV A[ 



A. The Star Graph Model 

The star graph is numerically modeled as a set of in- 
terconnected transmission lines as shown, schematically, 
in Fig. [9] [36]. This is a one port system, hence, the input 
signal is injected into the driving transmission line and 
the output signal is also retrieved from the same line. 
The driving transmission line has zero length. The driv- 
ing transmission line is connected with a number of trans- 
mission lines which are all connected in parallel with each 
other. The transmission line properties (for a line labeled 
by i) are length (Li), characteristic admittance (Y C i), fre- 
quency dependent complex propagation constant (7^(0;)), 
and complex reflection coefficient (r^) (for reflection from 
the terminations of the lines that are not connected to the 
driving line). The driving line has zero length, thus its 
only adjustable property is its characteristic admittance 
(Y cd ). ' 

This one port system is modeled by using the analyti- 
cally derived expression for its scattering parameter as a 
function of frequency, Sn(u). The scattering parameter 
can be expressed in terms of the characteristic admit- 
tance of the driving line (Y ca \) and the input admittance 
(Yi) of each of the other transmission lines while looking 
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FIG. 9: Schematic of the star graph model. There are N 
transmission lines that are connected in parallel, and a driving 
transmission line of zero length. Each of the N lines (labeled 
by i) can have unique length (Li), characteristic admittance 
(Y C i), frequency dependent propagation constant (7^(0;)), and 
reflection coefficient (1^). The driving line has a characteristic 
admittance of Y c d . 



towards them, 



Ycd + EiYiiu) 



Su(u) 



(2) 



where N is the number of transmission lines connected 
in parallel. The input admittance of each transmission 
line (labeled by z), Yi{uo), can be expressed in terms of 
the above mentioned properties of the line, 



1 _ F . p 2-fi{u:)Li 



(3) 



Once the scattering parameter of the system is com- 
puted over a broad frequency range, the response to 
any time domain input signal can be calculated. This 
is done by Fourier transforming the input signal to the 
frequency domain, multiplying it by the scattering pa- 
rameter (Sh(uj)) and inverse Fourier transforming the 
product back to the time domain. This establishes the 
star graph model as a time domain simulation of this 
quasi- ID system. However, in this section we are mainly 
interested in the frequency domain representation of the 
star graph using its scattering parameter. The frequency 
domain approach will be shown to be effective in quan- 
tifying strong perturbations that could not have been 
measured otherwise. 



B. Quantifying Strong Perturbations to the Star 
Graph 

As discussed in Sec. [TTJ the scattering parameter of a 
cavity can be used to quantify a VCSPP. The star graph 



is a quasi- ID system. A perturbation that changes the 
length of all of the lines of the star graph by the same 
proportion (i.e. with constant P) changes its effective 
volume while leaving its shape intact (i.e. it is a VCSPP). 
The effective volume of a star graph with N lines is given 
by 2 YZZT L i i 37 • Therefore, for a VCSPP perturbation 
that scales the lengths by a factor of P, the effective 
volume of the star graph also changes by a factor of P. 
This is diff erent from the case of the 3D cavity discussed 
in Sec. Ill (i.e. the volume of the 3D cavity changes by 
P 3 ). 

The baseline star graph was set up using the follow- 
ing parameter values. There were 10 lines whose length 
is given by Li = 0.3\10Vi] m (for i ranging from 1 to 
10). Each of these lines had a characteristic admittance, 
Yd, of IS. The characteristic admittance of the driv- 
ing line was chosen such that Y c d = X^=i° Yd = 105 in 
order to eliminate prompt reflection of signals injected 
through the driving line. The propagation constant of 
the lines is a function of the frequency, cj, and was given 
by 7i(co>) = 1^ with i — l,...,iV, where c is the speed 
of light in vacuum; thus, the lines themselves were con- 
sidered to be lossless (i.e. ) does not have a real 
part). However, energy was dissipated during reflection 
from the terminations of the lines. The reflection coef- 
ficient was given by I\ = e ct with i = l,...,iV, 
where r = 1.5/is. The amount of dissipation was de- 
signed to be independent of the size of the star graph. 
However, the dissipation introduced through sub-unitary 
values of I\ can be interpreted as an equivalent loss that 
could be introduced through ji(oo) (i.e. by introducing 
-^p as the real part of ji(uj) = 1^ + -^p). Thus, if the 
dissipation were modeled using 7i(cj), r = 1.5 {is would 
be interpreted as the time it takes the signals to decay 
by 1/e as they propagate along the lines. As a result, the 
typical 1/e decay time of the sona from the star graph 
was w 1.5/is, a value typical of our 3D experiment. 

The perturbed star graph was set up using identical 
values of parameters as the baseline star graph, except 
for the length, L^. The lengths of the perturbed star 
graph were chosen to be PL^, where P is the perturbation 
strength. The driving line has zero length in both the 
baseline and perturbed star graphs. 



1. Using the Time Domain Approaches 

A Gaussian pulse of width Ins and center frequency 
1GHz was used to generate a baseline and perturbed 
sona from the baseline and perturbed star graphs. The 
time domain sensing technique which is based on scat- 
tering fidelity was applied to quantify a perturbation of 
strength P = 1.0002. Fig. [10| shows the scattering fidelity 
of the baseline and perturbed sonas before (blue) and af- 
ter (red) optimally stretching the baseline sona. This 
demonstrates that the scattering fidelity sensing tech- 
nique can be used to quantify a VCSPP in the quasi- ID 
chaotic system. The result gives the clearest evidence 
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FIG. 10: Scattering fidelity of baseline and perturbed sonas 
from the star graph before (blue .) and after (red +) optimally 
scaling the baseline sona along its time axis. The perturbation 
is a VCSPP with P = 1.0002. 




FIG. 12: The average SF(t) of the baseline sona and per- 
turbed sonas, which are optimally scaled by P (perturbation 
magnitude) along their time axis. The sonas were collected 
from the baseline star graph, and a perturbed star graph (VC- 
SPP with scaling factor P). The average SF(t) was taken 
from time t = 0s to t = 10/is of the sonas; error bars show 
the associated standard deviation. 




1.0001 1.0002 1.0003 

Scaling Factor 



1.0004 



FIG. 11: Peak-to-peak amplitude (PPA) of time reversed 
pulse reconstructed inside the perturbed star graph using the 
baseline sona which is scaled along its time axis by different 
scaling factors. 



to the discussi on of fidelity decay induced by VCSPPs 
in Sec. IIID1 (i.e. the results shown in Figs. [5] and [6]). 
VCSPPs induce a scattering fidelity oscillation that is 
superimposed on a fidelity decay (see Fig. [l0|). 

The application of the time domain sensing technique 
which is based on time reversal is presented in Fig. [TT] 



Fig. 11 shows the PPA of the time reversed pulse recon- 
structed using the baseline sona scaled with different fac- 
tors along its time axis. The optimal PPA was obtained 
when the sona was scaled by a factor exactly equal to the 
stretching of the transmission line lengths (i.e. 1.0002). 

The perturbation strength of P = 1.0002 was shown to 
be detectable using time domain techniques in the star 



graph set up described. However, as the perturbation got 
stronger, the shortcoming of the time domain techniques 
was revealed. For perturbation strength values, P, rang- 
ing from 1.0002 to 1.14, the SF(t) of the perturbed sona 
and the optimally scaled (along the time axis) baseline 
sona was examined. For each P, the SF(t) was averaged 
over time, t (from Os to 10/is, which is the duration of the 
sonas). The closeness of the average SF(t) to 1 indicates 
the success of the time domain technique to undo the 
effect of the VCSPP, and to measure it. Fig. [12] shows 
the average value of SF(t) (for optimally scaled baseline 
sona) versus log 10 (P); the standard deviations of SF(t) 
taken over t are also shown as an error bar in Fig. [l2| 
As P increased beyond about 10 025 « 1.06, the effec- 
tiveness of the time domain sensing technique starts to 
deteriorate. 



Based on the discussion in Sec. |III E[ VCSPP pertur- 
bations that are characterized by P > ^ max ~ q'Ighz ~ 
1.15 are not expected to be quantified using the time 
domain techniques. The mean spacing between reso- 
nant frequencies of the baseline star graph is Af = 
^iZ N - « 2MHz, where N is the number of lines in 

the star graph [37], and f max - f min = 1GHz > Af. 
Therefore, the probing pulse excites several resonances 
of the baseline cavity. The result illustrated in Fig. [12 
shows that VCSPP perturbations that are characterized 
by P < 1.15 = 10 06 are quantifiable. 
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FIG. 13: Quantifying a volume changing perturbation in the 
frequency domain, (a) The |Sii| 2 (cj) of the baseline (blue) 
and a perturbed (red) star graph for perturbation strength 
P = 1.14. (b) The \Sh\ 2 (uj) of the baseline star graph (blue), 
and the |Sii| 2 (cj) of the perturbed star graph (black) after 
optimal frequency scaling which measures the VCSPP. 



V. DISCUSSION 

The development of this quantitative sensor opens up 
possibilities for several potentially useful applications. 
For example, when a cavity that has a homogenous mate- 
rial make up is cooled down (or warmed up) , it is interest- 
ing to check if the temperature stays uniform throughout 
all parts of the cavity. The sensor developed in this pa- 
per would allow one to check if the volume of the cavity 
is decreasing (or increasing) while the shape is intact; 
the sensor would also allow one to measure by how much 
the volume (and hence the temperature) of the cavity 
is changing. This is feasible as long as the temperature 
change primarily affects the volume of the cavity, and 
not the dielectric constant of the medium filling the cav- 
ity. This is one possible application of the VCSPP sen- 
sor which can practically compete with the traditional 
option of installing thermometers throughout the cav- 
ity. Another possible application of the VCSPP sensor 
is monitoring if a fluid has displaced another fluid uni- 
formly throughout a cavity. Since the speed of the waves 
inside different fluids can be different, the displacement 
of the fluid can change the volume of the cavity, as seen 
by the waves. This assumes that other wave properties 
(such as dissipation and dispersion) of the two fluids are 
similar. 



2. Using the Frequency Domain Approach 

VI. CONCLUSION 



The limitation of the time domain approach was dis- 
cussed in Sec. 



graph in Sec. I 



HIE 



VB1 



and demonstrated using the star 



As shown in Fig. 12 , for a strong 
perturbation such as P = 1.14, the time domain sens- 
ing techniques fail to undo the effect of the VCSPP, and 
hence to measure it. Here, a frequency domain approach 
is used to illustrate how the VCSPP with P = 1.14 can 
be measured. 



Fig. 13 'a) shows the |Sii| (oj) of the baseline (blue) 



and perturbed (red) star graphs. The |Sii| 2 (a;) of the 
perturbed star graph (larger in size by a factor of P = 
1.14) is compressed along the frequency axis compared to 
the baseline's case. This effect was predicted in Sec. [XTJ 
and schematically illustrated in Fig. [I] The frequency 
domain approach to measure VCSPP involves optimally 
scaling the frequency axis of the scattering parameter of 
the perturbed system (star graph in this case) to align 
it with the scattering parameter of the baseline system. 



Fig. Il3[b) shows the |5n| 2 (o;) of the baseline star graph 



(blue) and the optimally stretched |5ii| 2 (o;) of the per- 
turbed star graph (black) . The optimal frequency scaling 
factor was P = 1.14, which also successfully measures the 
VCSPP induced on the baseline star graph. To conclude, 
given the scattering parameters of a baseline and a per- 
turbed system, one can check if a VCSPP happened, and 
quantify it. 



Quantifying a generic perturbation to a wave chaotic 
enclosure is a challenge because the effect of a pertur- 
bation on the dynamics of waves is not generally very 
well understood. However, the effect of a perturbation 
that changes the volume but keeps the shape intact can 
be theoretically predicted. The theoretical prediction is 
most clear in the frequency domain. Thus, quantifying 
volume changing perturbations is best done in the fre- 
quency domain. Nonetheless, time domain approaches 
can be preferred for practical purposes. The time domain 
approach is limited by a maximum perturbation that can 
be measured. This limitation was demonstrated using a 
simulation of a star graph, which is a representative wave 
chaotic system. The time domain approach can work us- 
ing either scattering fidelity techniques or the classical 
analog of the Loschmidt echo. Quantification of a volume 
changing perturbation was experimentally demonstrated 
using these techniques. The volume changing perturba- 
tion was induced experimentally by changing the electri- 
cal volume of the cavity. The results of the experiment 
were compared with FDTD simulation results of the cav- 
ity and good agreement was found. 
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